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o 

We study a toy model for an interacting scalar field theory in which the fundamental excitations are confined in the sense 
1 of having unphysical, positivity-violating propagators, a fact tracing back to a decomposition of these in propagators with 

qq , complex conjugate mass poles (the so-called i-particles). Similar two-point functions show up in certain approaches to 

■ gluon or quark propagators in Yang-Mills gauge theories. We investigate the spectrum of our model and show that 
suitable composite operators may be constructed having a well-defined Kallen-Lehmann spectral representation, thus 
allowing for a particle interpretation. These physical excitations would correspond to the "mesons" of the model, the 
latter being bound states of two unphysical i-particles. The meson mass is explicitly estimated from the pole emerging 
in a resummed class of diagrams. 

Q-i 

1 Introduction 

The analytic study of the spectrum predicted by confining theories still remains a major theoretical challenge. 
In particular, despite the intense efforts dedicated to investigating Yang-Mills theories and QCD using different 
| nonperturbative tools, a complete understanding of the mechanism of confinement is not yet available, let stand 

ly-j ■ alone a clean-cut theoretical derivation of the spectrum of massive hadrons and glueballs from (almost) massless 

gluons and quarks. Such efforts can be appreciated from e.g. tools like Dyson-Schwinger equations |[T] [2) or sum 

■ rules approaches J3]|4][5]. 

During the last decade there have been many interesting developments in the direction of describing the dynam- 
ics of confined degrees of freedom in Yang-Mills theories. Both analytic and lattice methods were broadly used to 
investigate the two-point functions of the elementary fields of quarks and gluons, where the Landau gauge-fixing 

■ has been extensively employed. Recently, very precise lattice studies brought the community to a consensus. The 
gluon propagator differs sensibly from its perturbative form of the Faddeev-Popov type: in the deep infrared regime 
it is characterized by a clear positivity violation [6, 7] and achieves a non-vanishing value at zero momentum in 
D = 4, 3 dimensions JSJ [9] [TOj [TT] [12] [T3j |T4j [T5l, while it vanishes at the origin in momentum space in D = 2 

En ED. 

From a theoretical point of view, the infrared behavior of the gluon propagator has been object of a quite intense 
debate over many years, requiring a common effort from several groups. Different nonperturbative techniques 
have been employed, such as: Dyson-Schwinger equations |[T7l l2l [T8l[T9ll20ll2Tll22l . the renormalization group 
equations [23 24], effective field theory framework (25, 26], the Gribov-Zwanziger approach 112711281 and its more 
recent refined version ll29ll30ll3Tll32ll33l . and other approaches 0411351 . These efforts have culminated in the so- 
called decoupling solution for the gluon propagator, which turns out to be in very good agreement with the lattice 
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data in D = 4, 3 dimensions. Moreover, in D = 2, these methods coalesce towards a scaling type solution which, 
unlike the decoupling case, vanishes at the origin in momentum space, being in agreement with the numerical data 
inD = 2 BSE). 

One of the analytic expressions f 14, 38 , 39l[l5) that is capable of precisely describing lattice data in the infrared 
domain in D = 4, 3 dimensions and that allows for an analytic treatment of the calculation of Green's functions 
is provided by the Refined Gribov-Zwanziger (RGZ) scenario. The Gribov-Zwanziger (GZ) action l27l l28l is 
associated with a formulation of the gauge path integral that takes into account the presence of Gribov copies. 
Its refined version l29l [30l 1331 accounts for the presence of dimension 2 condensates which are dynamically 
generated by the restriction to the Gribov region, necessary to deal with the Gribov issue. Interestingly, the RGZ 
analytic expressions for the gluon propagator are in good accordance with the most recent lattice simulations and 
can be naturally decomposed into two quasi-degrees of freedom with complex-conjugated squared masses, which 
have been called i -particles l40l . 

Despite these significant advances in the description of propagators of confined elementary degrees of freedom, 
the emergence of a physical spectrum out of the confined elementary excitations still remains an open question 
from the analytic perspective. With the aim of taking advantage of these developments in the description of gluon 
propagators to shed light into the spectrum problem, encouraging estimates have been found for the spectrum of 
the lightest glueball states using analytic tools available in the RGZ scenario 14111 . Moreover, the inclusion of quark 
degrees of freedom in the RGZ action seems feasible and interesting results concerning also chiral symmetry and 
the meson spectrum are currently being investigated. 

It should be noticed that a usefuQ analytical form of the propagators is crucial in the construction of a frame- 
work that can assess these physical questions that are usually not amenable to analytical tools. 

In this paper we investigate the effects of interactions in a toy model based on the idea of confinement as 
implemented via the RGZ propagator decomposed into i-particles. Our goal is to understand how these complex 
modes may combine themselves to generate physical propagating modes in the correlation functions of suitable 
composite operators (= bound states). Of course, tackling this problem directly in QCD or even in pure Yang- 
Mills is an involved task, due to the non-Abelian complex nature of the interactions. It is therefore instructive to 
approach the problem through a toy model, obtaining a qualitative description of the mechanism of the formation 
of bound states in theories with unphysical elementary excitations. 

Previous descriptions in this line have used the propagators with lattice-fitted parameters as nonperturbative 
objects that carried information about confinement and estimates of the lightest glueball masses without including 
interactions explicitly have been obtained. In the current work, we shall combine the nonperturbative propagators 
with interactions defined in a toy model so that it is possible to resum a whole class of diagrams contributing to the 
two-point function of composite operators. The result for the interacting correlator displays a well-defined Kallen- 
Lehmann spectral representation. Moreover, we show explicitly the appearance of physical poles, associated with 
bound states of i-particles. This picture provides therefore a concrete qualitative understanding of the dynami- 
cal formation of the physical spectrum of a confining theory, defined in terms unphysical, positivity-violating, 
elementary fields. 

This paper is organized as follows. In the next section the confining toy model of i-particles is presented, 
including in detail our modeling of the interactions and the definition of the physical composite operators. In 
section [3] we describe the nonperturbative computation of the two-point function of a suitable physical composite 
operator. Results for the spectrum of bound states in dimensions D = 2, 3 and 4 are discussed in section [4] The 
last section is then dedicated to conclusions and final remarks. 



2 The confining scalar field theory 

In this section we describe the model to be studied in the following. We start by defining the free sector of the theory 
highlighting its confining properties and their formulation in terms of i-particles. Then we discuss the introduction 
of interactions, which should be ultimately responsible for generating a physical spectrum of bound states. 



1 . In the sense that one can compute quantities with it in a controllable fashion. 
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2.1 The free confining theory 

The model we will study has its free sector defined in such a way as to display a confining propagator of the type 
found in the RGZ framework I 



_ k 2 + m 2 R+ R_ 

D{ - ' = k 4 + k 2 (M 2 + m 2 ) + M 2 m 2 + 20 4 = k 2 + M 2 + k 2 + M 2 (2 ' 1} 

where M and m are real masses, 9 is a real massive parameter and 



M 2 = _ > ± ^ — . i _ 2g)4 (2.2) 



(M 2 + m 2 ) , l(M 

m 2 -M r 



where R± = ± ^ . The right-hand side of eg. (12. Il l suggests that -D(fc) encodes the propagation of two 

quasi-modes with potentially complex conjugated poles, if (M 2 ~ m 2 ) 2 < 86* 4 . This clearly unphysical feature is 
interpreted as a manifestation of confinement. 

The simplest model leading to this type of propagator is of the form 

S = /^(-^M» +2 J^)^ (2.3) 

from which one easily checks that the scalar field tp has a free propagator (^{k)il)(—k)) = D(k). This model can 
be rewritten in a local form through the introduction of auxiliary fields 

d D x + M 2 )ip - 4>{~d 2 + m 2 )4> + 6> 2 V# + 4>) + uj{-d 2 + m 2 )^ (2.4) 

where (cj>, </>) is a pair of bosonic complex conjugated fields and (ui, ui) is a pair of anticommuting fields. 

A version of this system was studied in 11401 with M = m = 0, corresponding to the original Gribov propagator. 
Since this is a quadratic action it can be cast in a complete diagonal form through a change of variables (tp, <fi, <£>) — > 

(v,x,v) 

S = J d D x ( \K-d 2 + + \v(-d 2 + M-)v - \ v (-d 2 + ™ 2 ) v + ^{-d 2 + m 2 )uj \ . (2.5) 

where A, 77 and V are real fields, with V being the imaginary part of </>. 

In order to simplify the analysis, in what follows we shall set M — m. In this case the fields U, ui and uj 
decouple and the action ( I2.5I ) reduces to 



8 = ]rz ( 7}H~d 2 + m 2 + iV26 2 )\ + ^v(-d 2 + rn 2 - iV29 2 ) V . (2.6) 

From this expression one immediately sees that the fields A and 77 correspond to the propagation of unphysical 
modes with complex masses m 2 ± i^/29 2 . These are the i-particles of the model, namely 

(A(fc)A(-fc)) - ']^~TZ~/^ r 2 (2 ' 7) 

Expression ( I2.6I ) describes a theory in which the fundamental excitations are not part of the physical spectrum. 
More precisely, it is not possible to analytically continue the propagators ( I2.8b to Minkowski space-time in order 
to obtain a well defined particle interpretation for the fundamental fields A and 77. In this sense, we might say that 
the model displays tree-level confinement. 



Note that, even though the 7-particles action (12. 6b has imaginary mass terms, it is Hermitian if we observe 
that A^ = 77, which follows from (jy = cf>. As the actions ( 12.51 ) and ( 12.3b are equivalent and the latter is clearly 
Hermitian, so should the former be. 
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2.2 Interactions 



Many properties of the action ( 12.6b were studied in BUI . In the present work we analyze the system including an 
interaction between the i-particles. Despite of the use of unphysical elementary modes, we shall be able to show 
that the interacting system does display a spectrum containing meson-like physical states. 

In principle there are many ways of introducing interactions in a model with the tree level propagator in eq. fl2.ll >. 



To that purpose, we first note that the action ( 12.3b is invariant under the replacement 2 — > — 2 . In the i- 
particles formulation ( 12.6b this symmetry can be understood by noting that in a path integral formulation A and rj 
are just dummy integration variables. In fact, the partition function is easily seen to be left invariant by 6 2 —> —9 2 
and A •(->■ 1). Therefore, we shall require that the interactions preserve this feature. 

The simplest interaction we can add to a scalar model is a quartic coupling. Thus we propose the following 
i-particles interacting model: 

S = J d D x y \K-d 2 + m 2 + iV26 2 )\ + \]{-d 2 + m 2 - iV29 2 ) V + gi(Ary) 2 + g 2 {\ 4 + ry 4 ) ^ , (2.9) 

where g\ and g 2 are real couplings. 

It is interesting to see how these interactions reflect on the original formulation ( 12.4b (with M = m). If we 
write (j> = -ks(U + iV) with U and V real fields, one realizes that V decouples from ip. In fact, it turns out that the 
sector of the original scalar model corresponding to (12. 9b is given by 

S = J d D x(^i>(-d 2 + m 2 )^-^U{-d 2 + m 2 )U + \f26 2 ij)U + gi(ipU) 2 + g 2 ^ + U 4 )) (2.10) 

where g\ = % — 3g 2 and g 2 = ^ + s k. The relation with the i-particles, A, r\, is provided by setting 

V=^=(X- V ). (2.11) 

Another way of defining the same physics is to start from the so-called replica model introduced in l42l . In this 
model we consider two copies of the same theory and couple them through a soft term that has the same effect 
as the imaginary masses discussed above. Consider for instance two scalar fields <f>i and (j> 2 , each described by a 
theory with quartic coupling. In this case, the replica model is given by the action 

S = y'd D xQ0 1 (-a 2 + m 2 )0 1 + i0 2 (-a 2 + m 2 )(/ )2 +5(^ + 02) + ^ 2 ^i02^) (2.12) 

The fact that a unique mass m as well as a unique quartic coupling g has been employed in expression ( 12.12b 
follows by demanding that the action displays the mirror symmetry l42l 



01 -> <t>2 , (t>2 -Kh. ■ (2.13) 

The last term in expression ( 12.12b . which contains the mass parameter 6 2 , implements a soft coupling between the 
two replica. In fact, this action can be cast into an i-particles formulation through the change of variables: 

02 = -L(A + t?) (2.14) 

giving 

d D x QA(-5 2 +m 2 + i\/26l 2 )A+i77(-a 2 + m 2 - t\/2e' 2 )?y + 3 5 (A7/) 2 + |(A 4 + ry 4 )^ , (2.15) 

which describes the same physics as expression ( 12. 9b , provided one sets g\ = 6g 2 in order to implement the mirror 
symmetry. 
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2.3 Physical operators 



The i -particles formulation shows clearly that the excitations corresponding to the elementary fields (A, rf) are 
unphysical. Nevertheless, following the construction outlined in ll40l l42l . physical states can be introduced by 
constructing suitable composite operators out of the fields (A, rf) which exhibit desirable analyticity properties, as 
encoded in the Kallen-Lehmann spectral representation. 

Such composite operators are obtained by requiring that the fields (A, 77) enter pairwise, i.e. the operator con- 
tains as many fields of the type A as of the type 77. This will ensure that in the corresponding correlation function 
only complex conjugate pairs of i-particles will propagate in the Feynman diagrams, a property which provides a 
good analytic structure. In practice, these operators can be obtained by requiring that their correlation function is 
left invariant by interchanging A and rj 11401 l42l . 

In the present case, the simplest example of a local composite operator with the required physical properties 
is 0(x) = X(x)rj(x). In fact, at lowest order in the interactions (gi — g-2 — 0), it is known that this operator has 
good analytical properties. In ll40l it was shown that the correlation function (0(k)0(— k))o has a well defined 
Kallen-Lehmann spectral representation: 

F D {k 2 ) = (0(k)0(-k)) Q D = f -^Tn " 7= " 7— = 1^ E2 ^ldT , 

K ' v v; v Jnhu J (2ir) D {k-p) 2 +m 2 -iV26 2 p 2 +m 2 +iV26 2 Jr r + k 2 

(2.16) 

where the spectral functions pjj have the following expressions, see l40l for detail^]: 

Pd=2{t) = - — 1 === , (2.17) 

2n Vr 2 - 86» 4 - 4m 2 r 

Pd=s{t) = ^-^=, (2.18) 



1 /, 86» 4 4to 2 
(47) 



The threshold to is in all cases given by 

T Q = 2 (m 2 + \/m 4 + 26> 4 ) . (2.20) 

It is important to observe that all spectral functions pu are positive in the corresponding range of integration. This 
is a necessary condition in order to interpret (12.16b as the propagation of a physical mode in Minkowski spaced 

Notice also that, in D = 4, the integral ( 12.16b is divergent in the ultraviolet, meaning that the correlation 
function has to be properly renormalized. As it is customary when dealing with spectral representation, this can 
be implemented be employing a subtracted dispersion relation l44l l45l l46l . In particular, in our case, it will be 
sufficient to perform only one subtraction and consider the subtracted expression J-D=4(k' 2 ) — J-d=4(0)- 



3 Nonperturbative propagator of the composite operator and physical spectrum 

The next step in establishing a physical spectrum for this confining theory is to compute the correlation functions 
of the composite operator 0(x) = X(x)r)(x) for the interacting case. We shall see that the interactions between 
i-particles will generate a physical bound state, i.e. a pole on the negative real axis in the complex k 2 plane. 
This composite correlator featuring a real pole may then be properly rotated to Minkowski space fl4l l45l |46*1 . 
corresponding to a physical propagating state of the theory. 

Since a singularity is not generated at any fixed order of perturbation theory, one must work in a nonperturbative 
framework. Our approach will be inspired by the standard resummation of bubble diagrams, which is an exact result 
in 0(A f )-symmetric scalar theories with a large number of field components N — > 00 11471 . It should be noticed, 



2. A recent numerical treatment of such correlation functions has been presented in 1431 . 

3. It is useful to remind that, once the Kallen-Lehmann spectral representation for a given correlation function has been obtained, the rotation 
from Euclidean to Minkowski is well defined I44ll45ll46l . 
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however, that the confining toy model discussed in this work is intrinsically different from 0(N) scalar theories, 
since it involves necessarily a doublet of i-particles, presenting different (complex-conjugated) masses. 

The bubble diagram resummation for the composite correlation function (X(x)r)(x)X(y)r)(y)) in the confining 
theory defined in Eq.( l2.9l l is an exact result if one fixes gi = 0, ignoring quartic self interactions. Setting 92 = 
is in principle a well-defined procedure in dimensions D = 2 and 3, since the resulting theories are UV finite 
and stable under quantum corrections. In the D = 4 case, (divergent) quartic self interactions are generated via 
radiative corrections and these terms are in principle needed to guarantee renormalizability. The toy model in 
D = 4 with 92 = that we will solve exactly is therefore either a fine-tuned renormalizable field theory or an 
effective low-energy theory, protected in the UV by an energy cutoffl In any case, the procedure to be presented 
in what follows illustrates successfully how an interacting confining theory defined in terms of unphysical degrees 
of freedom may generate a physical spectrum of bound states in different space time dimensions. 

Let us now compute the exact two-point function of the composite operator O = Ar/. As discussed above, we 
work with 92 = so that the only interaction vertex is g\ (A77) 2 . The series of diagrams contributing to (O(x)O(y)) 
is: 

(0(x)0(y)) = <T> + O0> + O00> +•••> ( 3 - 21 ) 



where the first diagram corresponds to the free case discussed in the previous section and defined as Tnik 2 ) in 
momentum space, cf. eqs. ( 12.161 ). Here, solid/dashed lines represent A/77 propagators, full dots correspond to g\- 
vertices and the empty dots make explicit the absence of interactions at the extreme points x, y, corresponding to 
the insertions of the operators O(x), 0(y). We omitted all self-energy corrections to the bubble propagators under 
the working hypothesis that these are already good approximations to the full nonperturbative propagators^. 

Since the n-th term in this series has the general form 

n bubbles 



OO O = /(0 c * (a; " !;) ^( fc2 )[^( fc2 )] n "' ^ 

the full two-point function can be written exactly as the result of a geometric series with ratio giJ 7 £>(k 2 ). In 
momentum space, we have 

00 TT 1 1 2 \ 

(O(k)O(-k)) = T D (k 2 ) [gi-Mfc 2 )]" - 1 _ g ^ D ( fc 2) ■ ( 3 - 23 ) 



n=0 

2 _ A/f2 



where the difference in the denominator points out to the possible existence of a pole for a given k = —M 

The explicit conditions for the existence of a pole in the result for the correlator in Eq.( l3.23b can be found by 
expanding around k 2 = —M 2 : 

T D (-M 2 )+0(k 2 + M 2 ) 

(o{k)o{-k)} - - - 5i<Fr) (_ M 2) _ gi r D (-M 2 ) (k 2 + W] + o ([fc 2 + m 2 ] 2 ) ' (3 ' 24) 

Therefore, if the conditions 

(i) l-giF D {-M 2 ) = 0, (3.25) 



(ii) R D {M 2 ) = - > , i.e. F' D {-M 2 ) < (3.26) 



4. One may also see our result in D = 4 as that of the analysis of the full theory (with #2 ^ 0), but coming from a Bethe-Salpeter approach 
with a bubble approximation for the kernel. 

5. This mimics the observation that the tree level RGZ propagators can account well for the lattice gluon propagator for example [38 15]. 
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are satisfied, then the two-point correlator of the composite operator O assumes the form of a physical pole at M 2 , 
with a positive residue Rd(M 2 ): 

(0(k)0(-k)) ^ (3.27) 
which corresponds, of course, to a well-defined Kallen-Lehmann spectral representation. 



4 Results and discussion 

Given the results for the spectral representation of the correlator (0(k)0(— k)) in the free case, Eqs. d2.16b - d2.19l l, 
one may systematically solve the pole conditions in Eqs. d3.25l ) and d3.26l l in order to obtain the physical spec- 
trum of the confining theory for each set of values for the parameters (D, gi,m, 9). All results are shown in this 
section as dimensionless ratios, with the threshold for two-particle production, ^/tq, being the mass unit used for 
normalization. 

In what follows we analyze three types of massive theories, namely (i) nonconfining (to ^ 0, 9 = 0); (ii) GZ 
confining (m = 0, 6 ^ 0) and (iii) RGZ confining (m ^ 0, 9 ^ 0). Bound state solutions are found for all cases in 
all dimensions investigated, D = 2, 3 and 4. 




-M 2 

T Figure 2: Residue of the pole as a function of 

the (negative) squared mass of the bound state 
Figure 1: Inverse coupling as a function of the (negative) in dimension 2. 

squared mass of the bound state in dimension 2. 

The mass M of the bound state is of course a function of the mass parameters m and 9 and the coupling gi, 
which has positive mass dimension in D < 4. Its specific form depends moreover on the space time dimension 
D. Nevertheless, the qualitative picture is essentially the same for the three theories analyzed, nonconfining and 
confining of the GZ- and RGZ-types, as shown in Figures [T] [3] and [5] This is an interesting feature of the bound 
state spectrum in this toy model: its static properties have no particular trace of the type of composites inside the 
"mesonic" state; both confining and nonconfining low-lying excitations furnish qualitatively equivalent spectra, at 
least statically. 

In Figures [1] [3] and [5] one further verifies that the threshold for two-particle production, ^/tq, is a mass upper 
limit for the bound states found. The heaviest bound states are generated at very low coupling. As interactions 
are turned on, the binding energy increases and the bound state mass decreases accordingly. It is also shown in 
Figures [T] [3] and [5] that the existence of a bound state with mass which is considerably lower than the threshold 
^/tq requires large couplings. For instance, a bound state with half the mass threshold ^/tq requires couplings of 
(/i/tq ^ 2 ~ 2, 10 and 500 in dimensions D = 2,3 and 4, respectively. Furthermore, nearly massless bound 
states are found in the deeply nonperturbative regime. 

In addition, the results show that the theory is physically meaningful only for sufficiently small couplings, 
.9i < ffcrit, with g CI a = 1/J r (fc 2 = 0) being associated with the massless bound state solution. For couplings 
larger than this critical value, tachyonic solutions appear, signaling the fact that the toy model becomes ill-defined 
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Figure 3: Inverse coupling as a function of the (negative) 
squared mass of the bound state in dimension 3. Here all three 
theories give the same curve in this normalization scheme. 
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Figure 4: Residue of the pole as a function of 
the (negative) squared mass of the bound state 
in dimension 3. 
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Figure 5: Inverse coupling as a function of the (negative) 
squared mass of the bound state in dimension 4. 
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Figure 6: Residue of the pole as a function of 
the (negative) squared mass of the bound state 
in dimension 4. 



in this region of its parameter space. The emergence of the tachyon might signal the instability of the considered 
vacuum, akin to what happens in the Gross-Neveu model where a condensation of the relevant composite operator 
takes place accordingly ll48l . Discussion of this would however lead us too far beyond the scope of this paper. 
This tachyon solution appears in all dimensions investigated, even though it is not made explicit in Figure [5] 
In dimension D = 4 the critical coupling is directly related to the subtraction scale and assumes therefore a 
running form, <7 cr ;t(A), with A being the subtraction energy scale. In Figure [5] A is fixed, for definiteness, so that 
ffcrit(A) — > oo. In dimension D = 2, also the massless bound state, i.e. gi = g alt should be treated with some 
care, as is usual for massless particles in dimension D = 2 1 49 53] [54). 

It is also instructive to investigate the behavior of the residue Rd{M 2 ) of the bound state poles found, since it 
is related to the probability amplitude of finding such a state. Figures [2] [4] and [6] show the results for dimensions 
D = 2, 3 and 4, respectively. In line with the expectation that the bound state should disappear in the absence of 
interactions, the residue go to zero as one approaches g\ — > 0. For D = 2 and 3 the residue grows monotonically 
as the mass M of the bound state is decreased, remaining finite also in the massless limit. In contrast, for D = 4 
the residue goes to zero as M — > 0, indicating that the probability of finding this state vanishes. This feature is, 
however, highly dependent on our subtraction choice for Td=a - which fixes FD=i{k 2 = 0) = -so that it 
hardly represents a general physical result, but it is rather a quantity that deserves a renormalization group/scheme 
improved analysis. 



8 



5 Conclusions 



The mechanism that generates the physical spectrum in a confining theory is an outstanding theoretical problem. 
Inspired by the advances achieved in the last years in the description of propagators of confined elementary particles 
in Yang-Mills theories and QCD as well as by the encouraging estimates of glueball masses obtained in BP . we 
adopted the R(GZ) scenario of confinement in which fundamental excitations have unphysical, positivity-violating 
propagators described by the combination of modes with complex-conjugated masses, the so-called i -particles. As 
detailed in the introduction, these RGZ propagators are in good agreement with very precise lattice data for the 
two-point functions of gluons in the deep infrared. 

In this paper we constructed a scalar quasi-particle toy model of these i-particles with the aim of investigat- 
ing the role of interactions in the dynamical formation of the physical spectrum. Using resummation techniques, 
we have shown that the interacting two-point correlator of suitable composite operators in this confining theory 
displays a well-defined Kallen-Lehmann spectral representation with a physical pole appearing dynamically. This 
picture provides therefore a concrete qualitative understanding of the dynamical formation of the physical spec- 
trum of a confining theory of the RGZ or GZ types, defined in terms unphysical, positivity-violating, elementary 
fields with complex mass. 

We have analyzed the theory in space time dimensions D = 2, 3 and 4, showing results for the behavior of the 
pole mass M of the physical bound states generated and the corresponding residua as the parameters of the model 
are varied. Our findings for these static properties of the composite spectrum of nonconfining and confining theories 
of two types (GZ and RGZ) have proven to be qualitatively similar. It should be noticed that this analysis is done 
in the Euclidean theory and no statement can be made concerning nonstatic properties of the bound state solutions 
found. We expect e.g. that real-time properties, such as widths and decay rates, of the spectrum of confining i- 
particle theories may be strikingly different from that of its nonconfining analog, with real elementary masses. 

Bound state formation has been investigated previously in nonconfining scalar theories in different space time 
dimensions. The theories investigated in the literature are in general not the one we have analyzed here. In dimen- 
sion D = 2, the absence of bound states was proven for weakly coupled A0 4 theories 1501 . Our results do not 
contradict these findings because we are dealing with a different model, with two fields A and r) that interact as 
(Ary) 2 , so that quartic self interactions are strictly absent. The consideration of the toy model with a nonzero quartic 
coupling gi (A 4 + rf) would involve extra classes of diagrams that could possibly bring the residue R3 (cf. Figure 
|2]i faster to zero in the weakly coupled domain (in the vicinity of M 2 = tq), featuring then the disappearance of 
bound states as seen in llBTfl . In dimension D = 3, the existence of a nonzero mass gap for A0 4 theories BP was 
also obtained in the weakly coupled regime, which is already consistent with our findings for the toy model. Bound 
states in the large N limit of A(^0i) 2 theories were studied in e.g. l52l . whose results are consistent with our find- 
ings. It should be noticed however that, even though the large N limit also involves the resummation of bubble 
diagrams, the toy model adopted here and its results are intrinsically different from those due to the inequivalent 
interaction vertices. Therefore, we believe our findings are consistent with the literature for nonconfining theories, 
besides contributing nontrivially to the investigation of bound state formation in confining theories. 
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